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Jet shapes have been measured in inclusive jet production in proton-proton collisions at ^/s = 
7 TeV using 3 pb~^ of data recorded by the ATLAS experiment at the LHC. Jets are reconstructed 
using the anti-fet algorithm with transverse momentum 30 GeV < pr < 600 GeV and rapidity in 
the region \y\ < 2.8. The data are corrected for detector effects and compared to several leading- 
order QCD matrix elements plus parton shower Monte Carlo predictions, including different sets 
of parameters tuned to model fragmentation processes and underlying event contributions in the 
final state. The measured jets become narrower with increasing jet transverse momentum and the 
jet shapes present a moderate jet rapidity dependence. Within QCD, the data test a variety of 
perturbative and non-perturbative effects. In particular, the data show sensitivity to the details of 
the parton shower, fragmentation, and underlying event models in the Monte Carlo generators. For 
an appropriate choice of the parameters used in these models, the data are well described. 

PACS numbers: 13.85.Ni, 13.85.Qk, 14.65.Ha, 87.18.Sn 

I. INTRODUCTION 

The study of the jet shapes f]\ in proton-proton collisions provides information about the details of the parton- 
to-jet fragmentation process, leading to collimated Hows of particles in the final state. The internal structure of 
sufficiently energetic jets is mainly dictated by the emission of multiple gluons from the primary parton, calculable in 
perturbative QCD (pQCD) [2|]. The shape of the jet depends on the type of partons (quark or gluon) that give rise 
to jets in the final state Q, and is also sensitive to non-perturbative fragmentation effects and underlying event (UE) 
contributions from the interaction between proton remnants. A proper modeling of the soft contributions is crucial 
for the understanding of jet production in hadron-hadron collisions and for the comparison of the jet cross section 
measurements with pQCD theoretical predictions 0, [1]. In addition, jet shape related observables have been recently 
proposed Q to search for new physics in event topologies with highly boosted particles in the final state decaying 
into multiple jets of particles. 

Jet shape measurements have previously been performed in pp Q, e^p 0> s-^d e+e^ 9] collisions. In this paper, 
measurements of differential and integrated jet shapes in proton-proton collisions at ^/s — 7 TeV are presented for 
the first time. The study uses data collected by the ATLAS experiment corresponding to 3 pb^^ of total integrated 
luminosity. The measurements are corrected for detector effects and compared to several Monte Carlo (MC) predictions 
based on pQCD leading-order (LO) matrix elements plus parton showers, and including different phenomenological 
models to describe fragmentation processes and UE contributions. 

The paper is organised as follows. The detector is described in the next section. Section 3 discusses the simulations 
used in the measurements, while Section 4 and Section 5 provide details on jet reconstruction and event selection, 
respectively. Jet shape observables are defined in Section 6. The procedure used to correct the measurements for 
detector effects is explained in Section 7, and the study of systematic uncertainties is discussed in Section 8. The jet 
shape measurements are presented in Section 9. Finally, Section 10 is devoted to summary and conclusions. 



II. EXPERIMENTAL SETUP 

The ATLAS detector [ToJ covers nearly the entire solid angle around the collision point with layers of tracking 
detectors, calorimeters, and muon chambers. For the measurements presented in this paper, the tracking system and 
calorimeters are of particular importance. 

The ATLAS inner detector has full coverage in [ll| and covers the pseudorapidity range \ri\ < 2.5. It consists 
of a silicon pixel detector, a silicon microstrip detector and a transition radiation tracker, all immersed in a 2 Tesla 
magnetic field. High granularity liquid-argon (LAr) electromagnetic sampling calorimeters cover the pseudorapidity 
range \ri\ < 3.2. The hadronic calorimetry in the range I77I < 1.7 is provided by a scintillator-tile calorimeter, which is 
separated into a large barrel and two smaller extended barrel cylinders, one on either side of the central barrel. In the 
end-caps {\ri\ > 1.5), LAr hadronic calorimeters match the outer \r]\ limits of the end-cap electromagnetic calorimeters. 
The LAr forward calorimeters provide both electromagnetic and hadronic energy measurements, and they extend the 
coverage to |ry| < 4.9. 
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The trigger system uses three consecutive trigger levels to select events. The Level- 1 (LI) trigger is based on 
custom-built hardware to process the incoming data with a fixed latency of 2.5 us. This is the only trigger level 
used in this analysis. The events studied here are selected either by the system of minimum-bias trigger scintillators 
(MBTS) or by the calorimeter trigger. The MBTS detector fl^] consists of 32 scintillator counters of thickness 2 cm 
organized in two disks. The disks are installed on the inner face of the end-cap calorimeter cryostats ed z — ±356 cm, 
such that the disk surface is perpendicular to the beam direction. This leads to a coverage of 2.09 < jryj < 3.84. The 
jet trigger is based on the selection of jets according to their transverse energy, Ex- The LI jet reconstruction uses 
the so called jet elements, which are made of electromagnetic and hadronic cells grouped together with a granularity 
of A(j) X Ary = 0.2 x 0.2 for \ri\ < 3.2. The jet finding is based on a sliding window algorithm with steps of one jet 
element, and the jet Et is computed in a window of configurable size around the jet. 



III. MONTE CARLO SIMULATION 



Monte Carlo simulated samples are used to determine and correct for detector effects, and to estimate part of the 
systematic uncertainties on the measured jet shapes. Samples of inclusive jet events in proton-proton collisions at 
= 7 TeV are produced using both PYTHIA 6.4.21 [l^ and HERWIG++ 2.4.2 [3 event generators. These MC 
programs implement LO pQCD matrix elements for 2 — > 2 processes plus parton shower in the leading logarithmic 
approximation, and the string [Tsl and cluster (l6| models for fragmentation into hadrons, respectively. In the case 
of PYTHIA, different MC samples with slightly different parton shower and UE modeling in the final state are 
considered. The samples are generated using three tuned sets of parameters denoted as ATLAS-MC09 [pTj , DW [l^ , 
and Perugia2010 [19]. In addition, a special PYTHIA- Perugia2010 sample without UE contributions is generated. 
Finally, inclusive jet samples are also produced using the ALPGEN 2.13 |20| event generator interfaced with HERWIG 
6.5 dill and JIMMY 3.41 ^ to model the UE contributions. HERWIG-h-h and PYTHIA-MC09 samples are generated 
with MRST2007LO* [H parton density functions (PDFs) inside the proton, PYTHIA-Perugia2010 and PYTHIA-DW 
with CTEQ5L [24] PDFs, and ALPGEN with CTEQ61L 25] PDFs. 

The MC generated samples are passed through a full simulation [2^ of the ATLAS detector and trigger, based 
on GEANT4 [27]. The Quark Gluon String Precompound (QGSP) model [llj is used for the fragmentation of the 
nucleus, and the Bertini cascade (BERT) model [231 for the description of the interactions of the hadrons in the 
medium of the nucleus. Test-beam measurements for single pions have shown that these simulation settings best 
describe the response and resolution in the barrel i30|] and end-cap !31j calorimeters. The simulated events are then 
reconstructed and analyzed with the same analysis chain as for the data, and the same trigger and event selection 
criteria. 



IV. JET RECONSTRUCTION 

Jets are defined using the anti-Zct jet algorithm [sl] with distance parameter (in y—(f> space) R = 0.6, and the energy 
depositions in calorimeter clusters as input in both data and MC events. Topological clusters ^ are built around 
seed calorimeter cells with |£'coii| > 4cr, where cr is defined as the RMS of the cell energy noise distribution, to which 
all directly neighboring cells are added. Further neighbors of neighbors are iteratively added for all cells with signals 
above a secondary threshold ji^ceiil > 2ct, and the clusters are set massless. In addition, in the simulated events jets 
are also defined at the particle level [33^ using as input all the final state particles from the MC generation. 

The anti-kt algorithm constructs, for each input object (either energy cluster or particle) i, the quantities dij and 
diB as follows: 

= min(fc,f )^^i^, (1) 
d^B = k-^. (2) 

where 

{AR)l^{y,~y,f + {^,~<f>,)\ (3) 

kti is the transverse momentum of object i with respect to the beam direction, (pi its azimuthal angle, and yi its 
rapidity. A list containing all the dij and diB values is compiled. If the smallest entry is a dij, objects i and j are 
combined (their four- vectors are added) and the list is updated. If the smallest entry is a diB, this object is considered 
a complete "jet" and is removed from the list. As defined above, dij is a distance measure between two objects, and 



3 



Trigger Information 


PT (GeV) 


trigger configurations 


integrated luminosity (nb ^) 


30 - 60 


MBTS 


0.7 


60 - 80 


L1_5/MBTS 


17 


80 - 110 


L1_10/L1_5/MBTS 


96 


110 - 160 


L1_15/L1_10/L1_5/MBTS 


545 


160 - 210 


L1_30/L1_15/L1_10/L1_5/MBTS 


1878 


210 - 600 


L1_55/L1_30/L1_15/L1_10/L1_5/MBTS 


2993 



TABLE I: For the various jet pr ranges, the trigger configurations used to collect the data and the corresponding total 
integrated luminosity. MBTS denotes the use of the minimum-bias trigger scintillators, while Ll_5, Ll_10, Ll_15, Ll_30, and 
Ll_55 correspond to LI calorimeter triggers with 5, 10, 15, 30, and 55 GeV thresholds, respectively. 



diB is a similar distance between the object and the beam. Thus the variable i? is a resolution parameter which sets 
the relative distance at which jets are resolved from each other as compared to the beam. The anti-fct algorithm is 
theoretically well-motivated 32] and produces geometrically well-defined ( "cone- like" ) jets. 

According to MC simulation, the measured jet angular variables, y and 0, are reconstructed with a resolution 
of better than 0.05 units, which iinproves as the jet transverse momentum, px, increases. The measured jet px is 
corrected to the particle level scale [5| using an average correction, computed as a function of jet transverse momentum 
and pseudorapidity, and extracted from MC simulation. 



V. EVENT SELECTION 



The data were collected during the first LHC run at ^/s = 7 TeV with the ATLAS tracking detectors, calorimeters 
and magnets operating at nominal conditions. Events are selected online using different LI trigger configurations 
in such a way that, in the kinematic range for the jets considered in this study (see below), the trigger selection is 
fully efficient and does not introduce any significant bias in the measured jet shapes. Table 1 presents the trigger 
configurations employed in each px region and the corresponding integrated luminosity. The unprescaled trigger 
thresholds were increased with time to keep pace with the LHC instantaneous luminosity evolution. For jet pt 
smaller than 60 GeV, the data are selected using the signals from the MBTS detectors on either side of the interaction 
point. Only events in which the MBTS recorded one or more counters above threshold on at least one side are 
retained. For larger pr, the events are selected using either MBTS or LI calorimeter based triggers (see Section 2) 
with a minimum transverse energy threshold at the electromagnetic scale [3j| that varies between 5 GeV (Ll_5) and 
55 GeV (Ll_55), depending on when the data were collected and the pr range considered (see Table 1). 

The events are required to have one and only one reconstructed primary vertex with a z-position within 10 cm of 
the origin of the coordinate system, which suppresses pile-up contributions from multiple proton-proton interactions 
in the same bunch crossing, beam-related backgrounds and cosmic rays. In this analysis, events are required to 
have at least one jet with corrected transverse momentum pT > 30 GeV and rapidity \y\ < 2.8. This corresponds 
approximately to the kinematic region, in the absolute four momentum transfer squared - Bjorken-x plane, of 
10'^ GeV^ < < 4 X 10^ GeV^ and 6 x 10~^ < a; < 2 x 10~^. Additional quality criteria are applied to ensure that 
jets are not produced by noisy calorimeter cells, and to avoid problematic detector regions. 

VI. JET SHAPE DEFINITION 

The internal structure of the jet is studied in terms of the differential and integrated jet shapes, as reconstructed 
using the uncorrected energy clusters in the calorimeter associated with the jet. The differential jet shape p(r) as 
a function of the distance r = a/ Ay^ -|- A(/)^ to the jet axis is defined as the average fraction of the jet px that lies 
inside an annulus of inner radius r — Ar/2 and outer radius r + Ar/2 around the jet axis: 

11 ^prir - Ar/2,r + Ar/2) , „ . ,„ 

Pi^) - A-rl^Il f' ^[n) ' A./2<.<i?-A./2, (4) 

where pt(''i, ^2) denotes the summed pt of the clusters in the annulus between radius ri and r2, N^°^ is the number of 
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jets, and R = 0.6 and Ar = 0.1 are used. The points from the differential jet shape at different r values are correlated 
since, by definition, p{r) Ar — 1. Alternatively, the integrated jet shape 5'(r) is defined as the average fraction 
of the jet pt that lies inside a cone of radius r concentric with the jet cone: 



where, by definition, '^{r — R) = 1, and the points at different r values are correlated. The same definitions apply to 
simulated calorimeter clusters and final-state particles in the MC generated events to define differential and integrated 
jet shapes at the calorimeter and particle levels, respectively. The jet shape measurements are performed in different 
regions of jet px and |y|, and a minimum of 100 jets in data are required in each region to limit the statistical 
fiuctuations on the measured values. 



The measured differential and integrated jet shapes, as determined by using calorimeter topological clusters, are 
corrected for detector effects back to the particle level. This is done using MC simulated events and a bin-by-bin 
correction procedure that also accounts for the efficiency of the selection criteria and of the jet reconstruction in 
the calorimeter. PYTHIA-Perugia2010 provides a reasonable description of the measured jet shapes in all regions of 
jet Pt and |y|, and is therefore used to compute the correction factors. Here, the method is described in detail for 
the differential case. A similar procedure is employed to correct independently the integrated measurements. The 
correction factors U{r,pT, \ y\) are computed separately in each jet px and \y\ region. They are defined as the ratio 
between the jet shapes at the particle level p{r)P^J^, obtained using particle-level jets in the kinematic range under 
consideration, and the reconstructed jet shapes at the calorimeter level p{r)^j., after the selection criteria are applied 
and using calorimeter-level jets in the given pT and \y\ range. The correction factors U{r,pT, \y\) = p{r)^c / P{''')'mc 
present a moderate p^ and \y\ dependence and vary between 0.95 and 1.1 as r increases. For the integrated jet shapes, 
the correction factors differ from unity by less than 5%. The corrected jet shape measurements in each pT and \y\ 
region are computed by multiplying bin-by-bin the measured uncorrected jet shapes in data by the corresponding 
correction factors. 



A detailed study of systematic uncertainties on the measured differential and integrated jet shapes has been per- 
formed. The impact on the differential measurements is described here in detail. 

• The absolute energy scale of the individual clusters belonging to the jet is varied in the data according to studies 
using isolated tracks [5], which parametrize the uncertainty on the calorimeter cluster energy as a function of 
Pt and rf of the cluster. This introduces a systematic uncertainty on the measured differential jet shapes that 
varies between 3% to 15% as r increases and constitutes the dominant systematic uncertainty in this analysis. 

• The systematic uncertainty on the measured jet shapes arising from the details of the model used to sim- 
ulate calorimeter showers in the MC events is studied. A different simulated sample is considered, where 
the FRITIOF [H plus BERT showering model is employed instead of the QGSP plus BERT model. 
FRITOF-hBERT provides the second best description of the test-beam results ^ after QGSP-hBERT. This 
introduces an uncertainty on the measured differential jet shapes that varies between 1% to 4%, and is approx- 
imately independent of pt and \y\. 

• The measured jet pt is varied by 2% to 8%, depending on pT and \y\, to account for the remaining uncertainty 
on the absolute jet energy scale Q , after removing contributions already accounted for and related to the energy 
of the single clusters and the calorimeter shower modeling, as discussed above. This introduces an uncertainty 
of about 3% to 5% in the measured differential jet shapes. 

• The 14% uncertainty on the jet energy resolution Q translates into a smaller than 2% effect on the measured 
differential jet shapes. 




(5) 



VII. CORRECTION FOR DETECTOR EFFECTS 



VIII. SYSTEMATIC UNCERTAINTIES 
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• The correction factors are recomputed using HERWIGH — h, which implements different parton shower, frag- 
mentation and UE models than PYTHIA, and compared to PYTHIA-Perugia2010. In addition, the correction 
factors are also computed using ALPGEN and PYTHIA-DW for px < 110 GeV, where these MC samples pro- 
vide a reasonable description of the uncorrected shapes in the data. The results from HERWIG-I--I- encompass 
the variations obtained using all the above generators and are conservatively adopted in all px and \y\ ranges to 
compute systematic uncertainties on the differential jet shapes. These uncertainties increase between 2% and 
10% with increasing r. 

• An additional 1% uncertainty on the differential measurements is included to account for deviations from unity 
(non-closure) in the bin-by-bin correction procedure when applied to a statistically independent MC sample. 

• No significant dependence on instantaneous luminosity is observed in the measured jet shapes, indicating that 
residual pile-up contributions are negligible after selecting events with only one reconstructed primary vertex. 

• It was verified that the presence of small dead calorimeter regions in the data does not affect the measured jet 



The different systematic uncertainties are added in quadrature to the statistical uncertainty to obtain the final result. 
The total uncertainty for differential jet shapes decreases with increasing px and varies typically between 3% and 
10% (10% and 20%) at r = 0.05 (r = 0.55). The total uncertainty is dominated by the systematic uncertainty, except 
at very large pr where the measurements are still statistically limited. In the case of the integrated measurements, 
the total systematic uncertainty varies between 10% and 2% (4% and 1%) at r = 0.1 (r = 0.3) as px increases, and 
vanishes as r approaches the edge of the jet cone. 

Finally, the jet shape analysis is also performed using either tracks from the inner detector inside the jet cone, 
as reconstructed using topological clusters; or calorimeter towers of fixed size 0.1 x 0.1 {y — (j> space) instead of 
topological clusters as input to the jet reconstruction algorithm. For the former, the measurements are limited to jets 
with \y\ < 1.9, as dictated by the tracking coverage and the chosen size of the jet. After the data are corrected back 
to particle level, the results from these alternative analyses are consistent with the nominal results, with maximum 
deviations in the differential measurements of about 2% (5%) at r=0.05 (r=:0.55), well within the quoted systematic 
uncertainties. 



The measurements presented in this article refer to differential and integrated jet shapes, p{r) and ^(r), corrected 
at the particle level and obtained for anti-fct jets with distance parameter R — 0.6 in the region \y\ < 2.8 and 
30 GeV < px < 600 GeV. The measurements are presented in separate bins of px and \y\. Tabulated values of the 
results are available in the Appendix and in Ref. [36]. 




Figures [T] to [3] show the measured differential jet shapes as a function of r in different px ranges. The dominant 
peak at small r indicates that the majority of the jet momentum is concentrated close to the jet axis. At low px, 
more than 80% of the transverse momentum is contained within a cone of radius r = 0.3 around the jet direction. 
This fraction increases up to 95% at very high px, showing that jets become narrower as px increases. This is also 
observed in Fig.[31 where the measured 1 — ^'(0.3), the fraction of the jet transverse momentum outside a fixed radius 
r = 0.3, decreases as a function oi px- 

The data are compared to predictions from HERWIG++, ALPGEN, PYTHIA-Perugia2010, and PYTHIA-MC09 
in Fig. [U to Fig. \M.a); and to predictions from PYTHIA-DW and PYTHIA-Perugia2010 with and without UE con- 
tributions in Fig. HKb). The jet shapes predicted by PYTHIA- Perugia2010 provide a reasonable description of the 
data, while HERWIG++ predicts broader jets than the data at low and very high px- The PYTHIA-DW predictions 
are in between PYTHIA-Perugia2010 and HERWIG-f-t- at low px and produce jets which are slightly narrower at 
high Px- ALPGEN is similar to PYTHIA- Perugia2010 at low px, but produces jets significantly narrower than the 
data at high px- PYTHIA-MG09 tends to produce narrower jets than the data in the whole kinematic range under 
study. The latter may be attributed to an inadequate modeling of the soft gluon radiation and UE contributions in 
PYTHIA-MC09 samples, in agreement with previous observations of the particle fiow activity in the final state p^ . 
Finally, Fig. Hljb) shows that PYTHIA-Perugia2010 without UE contributions predicts jets much narrower than the 
data at low px- This confirms the sensitivity of jet shape observables in the region px < 160 GeV to a proper 
description of the UE activity in the final state. 

The dependence on \y\ is shown in Fig.O where the measured jet shapes are presented separately in five different jet 
rapidity regions and different px bins, for jets with px < 400 GeV. At high px, the measured 1 — ^'(0.3) shape presents 
a mild \y\ dependence, indicating that the jets become slightly narrower in the forward regions. This tendency is 



shapes. 



IX. RESULTS 
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observed also in the various MC samples. Similarly, Figs. [6] and [7] present the measured 1 — 4' (0.3) as a function of 
Pt in the different \y\ regions compared to PYTHIA-Perugia2010 predictions. The result of tests to the data in 
Fig. [7] with respect to the predictions from the different MC generators are reported in Table 7, for each of the five 
rapidity regions. Here the different sources of systematic uncertainty are considered independent and fully correlated 
across pr bins (see Appendix). As already discussed, PYTHIA-Perugia2010 provides the best overall description of 
the data, while PYTHIA-Perugia2010 without UE contributions and ALPGEN show the largest discrepancies. 

Finally, and only for illustration, the typical shapes of quark- and gluon-initiated jets, as determined using events 
generated with PYTHIA-Perugia2010, are also shown in Figs. H] and [71 For this purpose, MC events are selected 
with at least two particle-level jets with pT > 30 GeV and \y\ < 2.8 in the final state. The two leading jets in this 
dijet sample are classified as quark-initiated or gluon-initiated jets by matching (in y — (j) space) their direction with 
one of the outgoing partons from the QCD 2 — >■ 2 hard process. At low pj- the measured jet shapes are similar to 
those from gluon-initiated jets, as expected from the dominance of hard processes with gluons in the final state. At 
high pt, where the impact of the UE contributions becomes smaller (see Fig. [4Kb)), the observed trend with px in the 
data is mainly attributed to a changing quark- and gluon-jet mixture in the final state, convoluted with perturbative 
QCD effects related to the running of the strong coupling. 

X. SUMMARY AND CONCLUSIONS 

In summary, jet shapes have been measured in inclusive jet production in proton-proton collisions at ^/s = 7 TeV 
using 3 pb""'^ of data recorded by the ATLAS experiment at the LHC. Jets are reconstructed using the anti-Zcj 
algorithm with distance parameter R = 0.6 in the kinematic region 30 GeV < pr < 600 GeV and |y| < 2.8. The data 
are corrected for detector effects and compared to different leading-order matrix elements plus parton shower MC 
predictions. The measured jets become narrower as the jet transverse momentum and rapidity increase, although with 
a rather mild rapidity dependence. The data are reasonably well described by PYTHIA-Perugia2010. HERWIG++ 
predicts jets slightly broader than the data, whereas ALPGEN interfaced with HERWIG and JIMMY, PYTHIA-DW, 
and PYTHIA-MC09 all predict jets narrower than the data. Within QCD, the data show sensitivity to a variety 
of perturbative and non-perturbative effects. The results reported in this paper indicate the potential of jet shape 
measurements at the LHC to constrain the current phenomenological models for soft giuon radiation, UE activity, 
and non-perturbative fragmentation processes in the final state. 
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FIG. 1: The measured differential jet shape, p{r), in inclusive jet production for jets with \y\ < 2.8 and 30 GeV < pt < 110 GeV 
is shown in different pr regions. Error bars indicate the statistical and systematic uncertainties added in quadrature. The 
predictions of PYTHIA-Perugia2010 (solid lines), HERWIG++ (dashed lines), ALPGEN interfaced with HERWIG and JIMMY 
(dotted lines), and PYTHIA-MC09 (dashed-dotted lines) are shown for comparison. 
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FIG. 2: The measured differential jet shape, p{r), in inclusive jet production for jets with \y\ < 2.8 and 110 GeV < pr < 310 GeV 
is shown in different pt regions. Error bars indicate the statistical and systematic uncertainties added in quadrature. The 
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Appendix A: Data Points and Correlation of Systematic Uncertainties 



Data for differential and integrated measurements are collected in Tables 2 to 6, which include a detailed description 
of the contributions from the different sources of systematic uncertainty, as discussed in Section 8. 

A test is performed to the data points in Tables 5 and 6 with respect to a given MC prediction, separately in 
each rapidity region. The systematic uncertainties are considered independent and fully correlated across pt bins, 
and the test is carried out according to the formula 



whore dj is the measured data point j, mcj{s) is the corresponding MC prediction, and ,s denotes the vector of standard 
deviations, Sj, for the different independent sources of systematic uncertainty. For each rapidity region considered, the 
sums above run over the total number of data points in pt and five independent sources of systematic uncertainty, and 
the is minimized with respect to s. Correlations among systematic uncertainties are taken into account in mcj{s). 
The results for the different MC predictions are collected in Table 7, and indicate that PYTHIA-Perugia2010 
provides the overall best description of the data. 
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± 0.004 


± 0.003 


± 0.001 


± 


0.004 


0.45 


0.243 


± 0.002 


± 0.026 




± 


0.025 




± 0.003 




± 0.005 


± 0.003 


^F 0.005 


± 


0.002 


0.55 


0.155 


± 0.001 


± 0.024 




± 


0.023 




± 0.002 




± 0.002 


± 0.002 


^F 0.007 


± 


0.002 



TABLE II: The measured difTerential jet shape, p(r), as a function of r in different pr regions, 
[2]). The contributions from the different sources of systematic uncertainty are listed separately. 



for jets with \y\ < 2.8 and 30 GeV < pr < 210 GeV (see Figs. [T] and 



pjr) (0 < \y\ < 2.8) 













210 GcV < p^, < 


260 GcV 












P ± 


(stat.) ± 


(syst.) 




cluster c-scalc 


shower model 


jet c-scalc 


resolution 


correction 


non-closure 


0.05 


7.004 ± 0.021 ± 0.146 




T 


0.061 


± 0.084 




T 0.037 


± 0.055 


± 0.035 


± 0.070 


0.15 


1.612 


± 0.012 


± 0.066 




T 


0.001 


T 0.050 




± 0.019 


± 0.033 


^ 0.012 


± 


0.016 


0.25 


0.672 


± 0.008 


± 0.035 




± 


0.012 


T 0.027 




± 0.011 


± 0.013 


T 0.005 


± 


0.007 


0.35 


0.353 


± 0.005 


± 0.024 




± 


0.019 


T 0.010 




± 0.007 


± 0.006 


:p 0.005 


± 


0.004 


0.45 


0.212 


± 0.003 


± 0.024 




± 


0.020 


T 0.007 




± 0.005 


± 0.004 


T 0.008 


± 


0.002 


0.55 


0.136 


± 0.001 


± 0.020 




± 


0.019 


T 0.001 




± 0.003 


± 0.003 


T 0.005 


± 


0.001 














260 GcV < p^, < 


310 GcV 












P ± 


(stat.) ± 


(syst.) 




cluster c-scalo 


shower model 


jet c-scalc 


resolution 


correction 


non-closure 


0.05 


7.300 


± 0.036 ± 0.113 




T 


0.053 


± 0.055 




T 0.027 


± 0.030 


T 0.001 


± 0.073 


0.15 


1.463 


± 0.021 


± 0.038 




± 


0.001 


T 0.030 




± 0.016 


± 0.009 


± 0.004 


± 


0.015 


0.25 


0.619 


± 0.014 


± 0.024 




± 


0.011 


T 0.013 




± 0.008 


± 0.013 


± 0.001 


± 


0.006 


0.35 


0.315 


± 0.008 


± 0.019 




± 


0.016 


T 0.010 




± 0.002 


± 0.003 


± 0.003 


± 


0.003 


0.45 


0.1S6 


± 0.004 


± 0.018 




± 


0.017 


T 0.005 




± 0.003 


± 0.002 


± 0.002 


± 


0.002 


0.55 


0.115 


± 0.002 


± 0.016 




± 


0.015 


^ 0.001 




± 0.002 


± 0.004 


T 0.004 


± 


0.001 














310 GcV < 
















P± 


(stat.) ± 


(syst.) 




cluster c-scalo 


shower u u n li 1 




resolution 


correction 


non-closure 


0.05 


7.495 


± 0.052 


± 0.128 




T 


0.043 


± 0.059 




^ 0.034 


± 0.033 


± 0.056 


± 


0.075 


0.15 


1.405 


± 0.031 


± 0.070 




± 


0.001 


T 0.056 




± 0.019 


± 0.026 


T 0.024 


± 


0.014 


0.25 


0.536 


± 0.018 


± 0.023 




± 


0.008 


T 0.008 




± 0.006 


± 0.001 


T 0.019 


± 


0.005 


0.35 


0.2S5 


± 0.011 


± 0.016 




± 


0.013 


± 0.002 




± 0.006 


± 0.004 


T 0.005 


± 


0.003 


0.45 


0.173 


± 0.006 


± 0.016 




± 


0.014 


T 0.001 




± 0.003 


± 0.005 


T 0.004 


± 


0.002 


0.55 


0.101 


± 0.003 


± 0.013 




± 


0.012 


± 0.001 




± 0.002 


± 0.001 


T 0.003 


± 


0.001 














400 GcV < P7, < 


500 GcV 












P ± 


(stat.) ± 


(syst.) 




cluster c-scalc 


shower model 


jet c-scalc 


resolution 


correction 


non-closure 


0.05 


7.720 


± 0.114 


± 0.106 




T 


0.034 


± 0.043 




T 0.031 


± 0.011 


± 0.036 


± 


0.077 


0.15 


1.339 


± 0.075 


± 0.054 




± 


0.001 


T 0.047 




± 0.020 


± 0.010 


T 0.003 


± 


0.013 


0.25 


0.489 


± 0.039 


± 0.023 




± 


0.006 


± 0.001 




± 0.008 


± 0.005 


T 0.020 


± 


0.005 


0.35 


0.226 


± 0.019 


± 0.012 




± 


0.009 


T 0.003 




± 0.003 


± 0.003 


T 0.005 


± 


0.002 


0.45 


0.128 


± 0.009 


± 0.011 




± 


0.010 


± 0.001 




± 0.002 


± 0.002 


T 0.003 


± 


0.001 


0.55 


0.086 


± 0.006 


± 0.011 




± 


0.010 


± 0.001 




± 0.001 


± 0.001 


T 0.003 


± 


0.001 














500 GcV < P7, < 


600 GcV 












P ± 


(stat.) ± 


(syst.) 




cluster c-scalc 


shower model 


jet c-scalc 


resolution 


correction 


non-closure 


0.05 


7.638 


± 0.261 ± 0.093 




T 


0.026 


± 0.001 




T 0.022 


± 0.009 


± 0.040 


± 0.076 


0.15 


1.400 


± 0.168 


± 0.037 




± 


0.001 


T 0.011 




± 0.013 


± 0.003 


T 0.030 


± 


0.014 


0.25 


0.475 


± 0.074 


± 0.017 




± 


0.006 


± 0.008 




± 0.009 


± 0.010 


T 0.003 


± 


0.005 


0.35 


0.257 


± 0.054 


± 0.012 




± 


0.010 


T 0.004 




± 0.004 


± 0.001 


T 0.002 


± 


0.003 


0.45 


0.153 


± 0.037 


± 0.012 




± 


0.011 


± 0.002 




± 0.004 


± 0.001 


T 0.002 


± 


0.002 


0.55 


0.078 


± 0.015 


± 0.010 




± 


0.009 


± 0.004 




± 0.002 


± 0.001 


^ 0.003 


± 


0.001 



TABLE III: The measured differential jet shape, p(r), as a function of r in different pr regions, for jets witli 
contributions from the different sources of systematic uncertainty are listed separately. 



< 2.8 and 210 GeV < pr < 600 GeV (see Fig.[3|. The 



1 - ^(r = 0.37 



















(0 < Ivl 


< 2.8) 














PT 


GoV) 


1 - = 0.3) ± (stat.) ± (syst.) 




cluster c-Bcalc 


show 


cr model 


jot 


c-scalc 




olution 






30 


- 40 


0.2193 


± 


0.0006 


± 


0.0325 




± 


0.0212 


± 


0.0001 


± 


0.0105 


± 


0.0057 


± 


0.0216 




40 


- 60 


0.1733 


± 


0.0008 


± 


0.0221 




± 


0.0177 


± 


0.0006 


± 


0.0070 


± 


0.0041 


± 


0.0104 




60 


- 80 


0.1347 


± 


0.0004 


± 


0.0157 




± 


0.0138 


± 


0.0010 


± 


0.0035 


± 


0.0017 


± 


0.0064 




SO 


- 110 


0.1146 


± 


0.0003 


± 


0.0117 




± 


0.0109 


± 


0.0005 


± 


0.0025 


± 


0.0007 


± 


0.0033 




110 


- 160 


0.0942 


± 


0.0003 


± 


0.0092 




± 


0.0084 


± 


0.0001 


± 


0.0021 


± 


0.0006 


± 


0.0030 




160 


- 210 


0.0789 


± 


0.0004 


± 


0.0067 




± 


0.0063 


± 


0.0007 


± 


0.0010 


± 


0.0008 




0.0015 




210 


- 260 


0.0698 


± 


0.0006 


± 


0.0069 




± 


0.0051 


± 


0.0015 


± 


0.0013 




0.0011 


± 


0.0020 




260 


- 310 


0.0616 


± 


0.0010 


± 


0.0046 




± 


0.0042 


± 


0.0014 


± 


0.0006 


± 


0.0008 


± 


0.0003 




310 


- 400 


0.0656 


± 


0.0015 


± 


0.0041 




± 


0.0035 


± 


0.0001 


± 


0.0010 


± 


0.0007 


± 


0.0016 




400 


- 500 


0.0442 


± 


0.0024 


± 


0.0033 




± 


0.0028 


± 


0.0001 


± 


0.0006 


± 


0.0006 


± 


0.0016 




500 


- 600 


0.0479 


± 


0.0070 


± 


0.0026 




± 


0.0022 


± 


0.0002 


± 


0.0008 


± 


0.0001 


± 


0.0012 





TABLE IV: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of pr, for jets with \y\ < 2.8 and 30 GeV < pr < 600 GeV (see Fig. |4}. The contributions 
from the different sources of systematic uncertainty are listed separately. 



1 - ^(r = 0.37 









(0 < y 


< 0.3) 














"jO V ) 


^ = 3) ± (stat ) ± (syst ) 






'ot c scale 




olution 


oorrcotion 










-1- n nt A n 
It u. u iiy 


ZIZ u.uuou 


11 2^ 


it 0.0057 


± 0.0234 






«n 




It U.Uioo 




It U.UU / 4 


-|- 




± 


0.0192 




«n 


an 


u.ioyo It u.uuii It u.UKt/ 


-I- n n 1 n 
It U.UiUD 


-J- n nn^n 
It u . UU / U 


It u.uuoy 


it 


0.0017 


± 


0.0062 




80 


' 110 


n 1 onQ -1- n nnriQ -1- n m i n 
u.i^uo It u.uuuy It u.uiiu 


It U . U U o D 


~\~ 0.0035 


It U.WZZ 


± 


0.0007 


± 


0.0055 




110 


- 160 


0.0990 ± 0.0007 ± 0.0087 


± 0.0067 


± 0.0025 


± 0.0017 


± 


0.0005 


it 0.0047 




160 


- 210 


0.0831 ± 0.0010 ± 0.0074 


± n.0051 


± 0.0004 


± 0.0008 


± 


0.0008 


it 


0.0053 




210 


- 260 


0.0758 it 0.0015 it 0.0047 


± n.0042 


± 0.0017 


± 0.0008 


± 


0.0011 


it 


0.0006 




260 


- 310 


0.0639 it 0.0024 ± 0.0068 


± 0.0035 


± 0.0032 


± 0.0003 


± 


0.0008 


± 


0.0049 




310 


- 400 


0.0578 it 0.0031 it 0.0034 


± 0.0030 


± 0.0002 


± 0.0013 


it 


0.0007 


it 


0.0007 




400 


- 500 


0.0486 ± 0.0044 ± 0.0037 


± 0.0024 


± 0.0022 


± 0.0006 


± 


0.0006 


± 


0.0017 










(0.3 < \y 


< 0.8) 












PT 


GoV) 


l-«(r = 0.3)± (atat.) ± (ayut.) 


cluster c-acalc 




jot o-scalo 




olution 






30 


- 40 


0.2219 it 0.0012 it 0.0390 


± 0.0173 


± 0.0036 


± 0.0109 


± 


0.0057 


± 


0.0326 




40 


- 60 


0.1779 it 0.0017 it 0.0233 


± 0.0145 


± 0.0051 


± 0.0059 


± 


0.0041 


± 


0.0160 




60 


- 80 


0.1378 it 0.0008 it 0.0159 


± 0.0117 


± 0.0021 


± 0.0041 


± 


0.0017 


it 


0.0097 




80 


110 


0.1179 it 0.0007 it 0.0116 


± 0.0093 


± 0.0002 


± 0.0025 


± 


0.0007 


it 


0.0063 




110 


- 160 


0.0963 it 0.0006 it 0.0094 


± 0.0073 


± 0.0006 


± 0.0018 


± 


0.0005 


it 


0.0056 




160 


- 210 


0.0847 ± 0.0007 it 0.0061 


± 0.0055 


± 0.0017 


± 0.0017 


it 


0.0008 


it 


0.0011 




210 


- 260 


0.0718 it 0.0012 it 0.0067 


± 0.0045 


± 0.0023 


± 0.0016 


it 


0.0011 


it 


0.0039 




260 


- 310 


0.0631 it 0.0019 it 0.0042 


± 0.0038 


± 0.0009 


± 0.0008 


it 


0.0008 


it 


0.0010 




310 


- 400 


0.0623 ± 0.0030 ± 0.0042 


± 0.0031 


± 0.0016 


± 0.0011 


± 


0.0007 


± 


0.0019 




400 


- 500 


0.0384 ± 0.0033 ± 0.0042 


± 0.0025 


± 0.0005 


± 0.0007 


± 


0.0006 


± 


0.0031 










(0.8 < |y 


< 1.2) 












PT 


GoV) 


1 - = 0.3) ± (stat.) ± (syst.) 


cluster c-acalc 


shower model 


jot o-scalo 




olution 






30 


- 40 


0.2191 ± 0.0014 ± 0.0314 


± 0.0233 


± 0.0030 


± 0.0102 


± 


0.0057 


± 


0.0172 




40 


- 60 


0.1736 ± 0.0020 it 0.0247 


± 0.0192 


± 0.0030 


± 0.0090 


± 


0.0041 


± 


0.0116 




60 


- 80 


0.1347 it 0.0009 it 0.0173 


± 0.0151 


± 0.0013 


± 0.0052 


± 


0.0017 


± 


0.0063 




80 


110 


0.1161 it 0.0008 it 0.0133 


± 0.0118 


± 0.0001 


± 0.0034 


it 


0.0007 


it 


0.0051 




110 


- 160 


0.0975 it 0.0007 it 0.0105 


± 0.0092 


± 0.0001 


± 0.0024 


it 


0.0005 


it 


0.0043 




160 


- 210 


0.0817 it 0.0009 ± 0.0071 


± 0.0069 


± 0.0007 


± 0.0014 


it 


0.0008 


it 


0.0010 




210 


- 260 


0.0721 it 0.0016 it 0.0073 


± 0.0054 


± 0.0010 


± 0.0015 


it 


0.0011 


it 


0.0044 




260 


- 310 


0.0639 it 0.0022 it 0.0051 


± 0.0046 


± 0.0010 


± 0.0016 


± 


0.0008 


± 


0.0002 




310 


- 400 


0.0529 ± 0.0031 ± 0.0058 


± 0.0038 


± 0.0001 


± 0.0009 


± 


0.0007 


± 


0.0042 




400 


- 500 


0.0593 ± 0.0079 ± 0.0037 


± 0.0030 


± 0.0014 


± 0.0011 


± 


0.0006 


± 


0.0012 





TABLE V: The measured integrated jet shape, 1 — ^'(r = 0.3), as a function of pr, for jets with 30 GeV < pr < 500 GeV in different jet rapidity regions (see Fig. [7|. 
The contributions from the different sources of systematic uncertainty are hsted separately. 



1 - ^(r = 0.37 











(1.2 < Ih 


< 2.1) 














PT 


GoV) 


l-«(r = 0.3)± (atat.) ± (ayot.) 


dust 


cr c-scalc 


shower model 


jot 


c-scale 




olution 






30 


- 40 


0.2177 ± 0.0010 ± 0.0325 


± 


0.0263 


± 0.0017 


± 


0.0114 


± 


0.0057 


± 


0.0140 




40 


- 60 


0.1731 ± 0.0014 ± 0.0244 


± 


0.0217 


± 0.0014 


± 


0.0066 


± 


0.0041 


± 


0.0077 




60 


- 80 


0.1331 ± 0.0007 ± 0.0178 


± 


0.0168 


± 0.0001 


± 


0.0035 


± 


0.0017 


± 


0.0046 




80 


110 


0.1130 ± 0.0006 ± 0.0140 


± 


0.0133 


± 0.0029 


± 


0.0029 


± 


0.0007 


± 


0.0012 




110 


- 160 


0.0904 ± 0.0005 ± 0.0109 


± 


0.0103 


± 0.0010 


± 


0.0019 


± 


0.0005 


± 


0.0029 




160 


- 210 


0.0735 ± 0.0007 ± 0.0082 


± 


0.0077 


± 0.0011 


± 


0.0015 


± 


0.0008 




0.0019 




210 


- 260 


0.0646 ± 0.0011 ± 0.0066 


± 


0.0061 


± 0.0007 


± 


0.0014 


± 


0.0011 


± 


0.0014 




260 


- 310 


0.0573 ± 0.0021 ± 0.0053 


± 


0.0051 


± 0.0007 


± 


0.0011 


± 


0.0008 


± 


0.0002 




310 


- 400 


0.0495 ± 0.0026 ± 0.0045 


± 


0.0043 


± 0.0005 


± 


0.0008 


± 


0.0007 


± 


0.0009 




400 


- 500 


0.0335 ± 0.0033 ± 0.0037 


± 


0.0035 


± 0.0006 


± 


0.0007 


± 


0.0006 


± 


0.0006 












(2.1 < Ijj 


< 2.8) 














PT 


GoV) 


1 - = 0.3) ± (stat.) ± (syst.) 


cluster c-scalc 




jet 






olution 






30 


- 40 


0.2110 ± 0.0014 ± 0.0256 


± 


0.0209 


± 0.0094 


± 


0.0098 


± 


0.0057 


± 


0.0003 




40 


- 60 


0.1664 ± 0.0021 ± 0.0193 


± 


0.0169 


± 0.0048 


± 


0.0066 


± 


0.0042 


± 


0.0023 




60 


- 80 


0.1274 ± 0.0011 ± 0.0153 


± 


0.0126 


± 0.0062 


± 


0.0057 


± 


0.0017 


± 


0.0012 




SO 


110 


0.1048 ± 0.0009 ± 0.0110 


± 


0.0099 


± 0.0031 


± 


0.0033 


± 


0.0007 


± 0.0004 




110 


- 160 


0.0S30 it 0.0008 ± 0.0090 


± 


0.0076 


± 0.0026 


± 


0.0034 


± 


0.0005 


± 


0.0019 




160 


- 210 


0.0626 ± 0.0010 ± 0.0074 


± 


0.0058 


± 0.0030 


± 


0.0026 


± 


0.0008 


± 


0.0020 




210 


- 260 


0.0607 ± 0.0023 ± 0.0066 


± 


0.0048 


± 0.0018 


± 


0.0027 


± 


0.0011 


± 


0.0029 




260 


- 310 


0.0538 ± 0.0040 ± 0.0047 


± 


0.0040 


± 0.0022 


± 


0.0007 


± 


0.0009 


± 


0.0006 





TABLE VI: The measured integrated jet shape, 1 — *I'(r — 0.3), as a function of pr, for jets with 30 GeV < pr < 500 GeV in different jet rapidity regions (see Fig. [7|. 
The contributions from the different sources of systematic uncertainty are hsted separately. 



to 
o 



x'/d.o.f 





0< \y\ <0.3 


0.3 < \y\ < 0.8 


0.8 < \y\ < 1.2 


1.2 < \y\ <2.1 


2.1 < \y\ < 2.8 


degrees of freedom (d.o.f) 


10 


10 


10 


10 


8 


PYTHIA-Perugia2010 


0.6 


1.8 


2.4 


1.4 


1.4 


HERWIG++ 


2.2 


2.3 


3.1 


1.8 


4.0 


PYTHIA-MC09 


1.0 


2.5 


2.4 


1.5 


3.2 


PYTHIA-DW 


2.4 


3.4 


6.9 


4.0 


5.2 


ALPGEN 


3.8 


9.8 


7.4 


6.7 


6.0 


PYTHIA-Perugia2010 (no UE) 


4.2 


9.7 


4.9 


8.6 


4.8 



TABLE VII: Results of tests to the data in Fig. [7] with respect to the different MC predictions. As discussed in the text, the different sources of systemati' 
uncertainty are considered independent and fully correlated across pr bins. 
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